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Abstract 

The superintegrability of so called Tremblay-Turbiner-Winternitz (TTW) 
model has been conjectured on the basis of the fact that all its trajectories are 
closed. This conjecture has been proven using the method based on solving the 
partial differential equations for two functions having the same Poisson bracket 
with the Hamiltonian. 

In the present short paper we show that superintegrability of TTW model can 
be established by using well-known elegant techniques of analytical mechanics. 
Moreover, the resulting expression ( after an appropriate ordering ) can be gen- 
eralized to the quantum-mechanical case. 
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I Introduction 



Recently an interesting integrable two-dimensional model has attracted some attention 
, both on classical and quantum levels [1] [3] - Its Hamiltonian reads 
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where without loosing generality one can put k > 0. The coordinate space is defined 
by the inequalities 0<r<oo,0<(^<^ (actually one can consider other sectors 

If < f < 2fc )• ^ ne study of classical trajectories generated by the Hamiltonian (jTJ) 
has been performed in Ref. [2J. It has been shown there that all bounded trajectories 
are closed for all rational values of k. Moreover, the period of motion is T = ^ so it 
is the same for all (bounded ) trajectories. 

These findings strongly suggest that the Hamiltonian ([1]) is superintegrable. This 
has been shown to be the case by Kalnins et al [4]. The authors of Ref. [4] used 
their own method based basically on solving the partial differential equations for two 
functions having the same Poisson bracket with the Hamiltonian so that their difference 
provides an additional constant of motion. 

The aim of the present short note is to show that the superintegrability of the model 
can be established using standard techniques of analytical mechanics [5], [6]. We show 
that the superintegrability takes place iff k is rational and find the explicit form of 
additional integral. 

II Superintegrability of TTW model for rational k 

The integrability of the model defined by the Hamiltonian (JTJ) is obvious because it 
admits separation of variables. Two independent commuting integrals of motions are 
the Hamiltonian itself as well as 
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In oder to find out whether and when our system is superintegrable we construct 
first the action variables. This is done using standard methods [5], [6]. The invariant 
Arnold - Liouville tori are given by the equations 



1 , 2 ak 2 13k 2 
_(V -| 1 _ — 

2 r cos 2 (kip) sin 2 (kip) 



p 2 r + co 2 r 2 + -(p 2 + — + r ) = E (3) 



2 ak 2 /3k 2 A 

Pi + Tn \ + —T7T^ = A ( 4 ) 



cos 2 (kip) sin 2 (kip) 



2 



As explained in Ref. [2], for bounded trajectories the following inequalities hold : E 2 — 
Au 2 > 0, A > 0, (A + (/3 - a)k 2 ) 2 - Ak 2 (5A > 0. 

Integrating over two generators of homotopy group of Arnold-Liouville torus one 
finds (after some trivial shift of I2 variable) the action variables ( cf. also [7]). 

I - E ^ (5) 



which leads to 



H = E = 4w(h + kl 2 ) (7) 

Therefore, calling the corresponding angle variables one obtains the equations 
of motion 
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Now, for k irrational the trajectory covers densely the invariant torus. Conse- 
quently, it cannot result from the intersection of the torus with some hypersurface 
corresponding to the constant value of additional globally defined integral of motion. 
The system is integrable but not superintegrable. 

Assume now that 



k=™ (10) 
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with m,n natural with no common divisor. It follows from eqs. ()8]) and (Q that 
rmfji — mp2 is an integral of motion which is functionally independent on Its 
disadvantage is that it is not single-valued function on phase space. In order to get a 
single- valued function one has to follow the standard procedure describe, for example, 
in [5], i.e. to take some periodic function of mipi — nip2, say cos(m-^i — n-0 2 ) or 
sin(m-0i — n^)- in practice, it is more convenient to consider 
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and then to take real or imaginary part. 

The angle variables are computed according to the formula [6] ( see also [7] ). 
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where S has the same meaning as in Ref. [2] once E and A are expressed in terms 
of action variables. Now, ^ can be readily obtained using eqs. (JSJ),® and ( IT2|) and 
the results of Ref.[2j ( eqs. (15), (19) and (20) therein ). There is no need to write out 
explicitly the form of angle variables. One only notes that consists of one term while 
■02 is a sum of three terms. In the combination rm/ji — n^, mifj% cancels against the 
first term in m/V The remaining contributions are computed from eqs. (19) and (20) 
of Ref. [2]. The only subtle point is that one has to find exp knowing sin if; which 
a priori is not unique. However, we know that exp (iifj) is single-valued on the phase 
space, so we expect that cos tp is determined from sin ip as a single-valued function; in 



fact, it appears that the relevant square root J 1 — sin 2 i() can be taken explicitly. 

Finally, writing out in explicit form eq. ffTTj) . we find the following additional integral 
of motion (actually, we simplified final expression by multiplying it by an appropriate 
function of E and A) : 



C= y ^ Pr +i(E-^yj (VA Pv sm2k(p + i{{P-a)k 2 + Acos2k V )) n (13) 

One can directly check that C is an integral of motion by taking the time derivative 
of C and using the Hamiltonian equations of motion or by computing the Poisson 
bracket of C with the Hamiltonian ([T]). This is completely straightforward. Taking real 
or imaginary part of C one obtains real integral of motion functionally independent 
on A and E. Obviously, we produce in this way only one new independent integral. 
Moreover, it is also straightforward to see that both real and imaginary parts contain 
only either even or odd positive powers of \J~A. Therefore, multiplying by y/A, if 
necessary, we obtain integrals polynomial in momenta p r and p v . One can do even 
better. Putting A = in eq. ffTB"]) one gets 

C(A = 0) =i m+n E m ((P- a) k 2 ) n (14) 

Therefore, the real (imaginary) part of C consists of integer powers of A for m + n even 
(resp. odd). By subtracting C(A = 0) and dividing by A one obtains a polynomial 
in momenta p r , p v of degree 2(m + n — l).The above result can be compared with 
the findings of Ref. 0| . Our final form of the integral is a polynomial in the momenta 
of degree smaller by one. This is because we made one step more. Namely, as it 
is explained above, the value of C(A = 0) (which is a function of energy ) has been 
subtructed which allowed to extract and neglect the A factor ( because it is the integral 
of motion by itself ). 
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The advantage of the formula (}Tlj) is that, under an appropriate ordering, C can 
be converted into its quantum counterpart see the forthcoming paper [Sj. It reads 



» _ r -sin(2fcy) d f y/A cos (2k<p) y/A f y/A 
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where A and v are defined by: a = A(A — 1) and j3 = v[y — 1) and trivial numeric 
factor has been omitted. Again, the fact that C is a constant of the motion can be 
verified by explicit computation of the commutator [C,H]. Moreover, C obeys the 
correspondence principle because, as it has been mentioned above, it is obtained from 
its classical counterpart by an appropriate ( although by far nontrivial ) ordering. 
It will be shown in Ref.([8]) that C (together with its hermitean conjugate) can be 
used to construct the polynomial (in derivatives) integral of motion. 



Ill Conclusions 

To sumarize let us discuss the relation of action-angle variables method used above to 
the one proposed in Refs.|l] and [9]. The common idea of both methods is that one 
is looking for the functions on phase space which develop linearly in time; once they 
are found their appropriate linear combination is a constant of motion. In fact, for 
2D systems admitting separation of variables the method of Refs.jl] and [5] provides 
a superintegral (i.e. a functionally independent additional integral of motion beyond 
those following from Liouville itegrability) which can also be obtained directly within 
Hamilton- Jacobi approach. However, in general this integral is not well-defined glob- 
ally on phase - space. This becomes crucial in the case of bounded motion (like in 
TTW model) which generically is quasiperiodic and, therefore, the additional integrals 
obtained in the above way are not separating ones. It is just a very idea of action- 
angle variables to normalize the new momenta such that the canonically conjugated 
coordinates are angles with 2tt - periodicity. This allows to decide at once whether 
well-defined superintegrals exist globally and to construct them by quadratures. More- 
over, within this framework, it is clear that the additional integrals do not need to be 
the polynomials in momenta. 
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